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Abstract 

Despite the recent interest in "organic spintronics" , the dominant spin relaxation mechanism 
of electrons or holes in an organic compound semiconductor has not been conclusively identified. 
There have been sporadic suggestions that it might be hyperfine interaction caused by background 
nuclear spins, but no confirmatory evidence to support this has ever been presented. Here, we 
report the electric-field dependence of the spin diffusion length in an organic spin-valve structure 
consisting of an Alq^ spacer layer, and argue that this data, as well as available data on the 
temperature dependence of this length, contradict the notion that hyperfine interactions relax 
spin. Instead, they suggest that the Elliott- Yafet mechanism, arising from spin-orbit interaction, 
is more likely the dominant spin relaxing mechanism. 
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Spin relaxation in most solids is caused primarily by mechanisms associated with spin- 
orbit and contact hyperfine interactions. Since compound organic semiconductors are typ- 
ically made of light elements (hydrogen, oxygen and carbon), the spin-orbit interaction in 
them should be very weak since it is proportional to the fourth power of the atomic num- 
ber of the constituent elements. At the same time, contact hyperfine interaction (between 
electron and nuclear spins) should also be very weak - at least in 7r-conjugated organic 
molecules - because the 7r-electrons' wavefunctions are mainly p z orbitals that have nodes in 
the molecular plane [lj. As a result, organics tend to exhibit long spin relaxation times that 
could exceed those in inorganic semiconductors by several orders of magnitude at temper- 
atures well above that of liquid nitrogen. This has generated significant interest in organic 
spintronics [2 5] owing to the realization that organic semiconductors could very well emerge 
as the material of choice in many spintronic applications. 

Despite all this interest and research, a question of fundamental importance has remained 
unanswered: which of the two mechanisms - spin-orbit interaction or hyperfine interaction 
- is the dominant causative agent for spin relaxation in organics. Because spin-orbit inter- 
action is so much weaker in organics than in inorganics, the trend has been to conjecture 
(tacitly) that contact hyperfine interaction must be the dominant spin relaxation mecha- 
nism 5H7J. However, to our knowledge, no conclusive evidence has ever been presented to 
substantiate this belief. This remains an open question. 

In this paper, we show that a large body of experimental evidence does not support the 
notion that hyperfine interaction is the dominant spin relaxation mechanism in the most 
widely studied organics. Instead, it points to the Elliott- Yafet mechanism js|, arising from 
spin-orbit interaction, as being the more likely culprit. We reach this conclusion based on 
the reported temperature and electric-field dependences of the spin diffusion length - the 
latter reported here - which are not consistent with hyperfine interaction, but are consistent 
with the Elliott- Yafet mechanism being dominant. In the rest of this paper, we elucidate 
the arguments leading to this conclusion. 

At a temperature T and in an electric field E, the spin diffusion length L S (E, T) of a spin 
carrier in any solid is related to the spin relaxation time t s (E,T) and the carrier mobility 
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n(E,T) according to the relation [12|, |l3] 



Ls(E,T) = 1 (1) 



2kT 1 V \2kT J 1 kTfi(E,T)Ts(E,T) 

where k is the Boltzmann constant and e is the electronic charge. 

Since carriers in organics travel by hopping from site to site assisted by thermal excitation 
in an electric field [similar to Poole Frenkel conduction] [l^] the mobility is usually expressed 



as 



15| 



fi(E, T) = tank ( ^] exp ( MVE-W} 

where d(E,T) is the mean hopping distance, t (E,T) is the mean hopping time, A(T) is 
the activation energy for hopping and (3(T) is the field emission constant. 
In low electric fields y\E\ <\J ef j,(E,T)r (et) )' Equation (1) simplifies to 



where (see Equation (2)) 



' kT 

[Ls(E,T)} low E « /i(0,T)r s (0,T), (3) 



Substituting Equation (4) in Equation (3), we get 



[WiL.,=««(o,r)^gg«p(^). (5) 

The last equation is very instructive. The quantities d(0,T), r (0,T) and A(T) should 
not have strong temperature dependence at cryogenic temperatures. Therefore, Equation 
(5) tells us that the low-field and low-temperature spin diffusion length [L S (E, T)] low E must 
increase with increasing temperature, unless r s (0, T) decreases with increasing temperature. 
This makes sense intuitively, and we could have predicted it from Equation (3) directly. 
Since carriers travel by Brownian motion, increased temperature should result in increased 
spin diffusion length unless the spin relaxation time decreases with increasing temperature. 
Thus, if we ever observe [L S (E, T)] low E decreasing with increasing temperature (at cryogenic 
temperatures), then we must conclude that r s (0, T) also decreases with rising temperature. 
Consequently, the temperature dependence of the spin diffusion length at low temperatures 



is very revealing; it tells us how the spin relaxation time varies with temperature. In turn, 
that can allow us to decipher the major spin relaxing mechanism. 

So far, every experiment reported in the literature has found that in organics, 
L S (E, T)] low E decreases with increasing temperature in the cryogenic range, either slowly 
9|, or moderately [10|, or rapidly That then tells us that r s (0,T) must also decrease 
with rising temperature, slowly, moderately or rapidly. The rapid decrease cannot be con- 
sistent with hyperfme interaction because in that mechanism, spin relaxation is caused by 
the magnetic field of the nuclear spins, which, at best, can have weak temperature depen- 
dence. Indeed, theories based on spin diffusion in a disordered organic in the presence of the 
ryperfine magnetic field predict a weak temperature dependence of the spin diffusion length 
6j. Therefore, the observation in ref. 11], which showed a rapid decrease of [L s (E,T)] low E 
with increasing temperature, is clearly inconsistent with the notion that hyperfme interac- 
tion could have been the dominant spin relaxing mechanism in that organic. 

To probe this matter further and correctly identify the dominant spin relaxation mecha- 
nism, we can investigate the electric-field dependence of the spin relaxation length and spin 
relaxation time in an organic, which, to our knowledge, has never been attempted. Here, we 
report some data on the electric-field dependence of the spin diffusion length in an organic 
and infer the electric-field dependence of the spin relaxation time from that data. This sheds 



further lig 



it on the dominant spin relaxation mechanism. 



In ref. 16], we carried out experiments in organic spin valves to extract the spin diffusion 



length under varying electric fields, where the organic was tris(8-hydroxyquinolinolato alu- 
minum) or Alqs with a chemical formula of C27H18N3O3AI. The spin valves were nanowires 
of Alq 3 with cobalt and nickel contacts. We will assume that Pi and P2 are the spin polariza- 
tions at the Fermi energy in the injecting and detecting contacts, a\ and 0.2 are the effective 
spin injection and detection efficiencies at the two contacts, and L is the organic layer thick- 
ness. Schottky barriers form at both contacts because of the energy level alignment 
The picture of carrier transport in these spin valves presented in refs. jsl, [lo| is that carriers 
first tunnel through the Schottky barrier at the injecting contact with a spin polarization 
P\di, then drift and diffuse through the bulk of the organic with exponentially decaying spin 
polarization e~ L ^ Ls( - E ' T \ and finally tunnel through the second Schottky barrier to reach the 
detecting contact. Therefore, the spin valve magnetoresistance ratio AR/R will be given by 
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the modified Julliere formula (adapted from ref. BQ): 

AR 2P 1 a 1 P 2 a 2 e- L/Ls(£ ' T) 



(6) 



R 1 - P 1 a 1 P 2 a 2 e- L / L °( E > T ) ' 

In order to verify the transport picture in refs. , we had measured the current- 

voltage (I-V) characteristics of the spin valves at varying temperatures {9]. They were nearly 
temperature-independent and almost piecewise linear. The current increased quasi-linearly 
with a small slope up to a threshold voltage of ~ 2 V, and then increased rapidly with 
a much larger slope. This behavior is inconsistent with tunneling through pinholes in the 
organic - which has been proposed as an alternate transport model jisj] - since that would 
have produced two features which are absent. First, tunneling causes the I-V characteristic 



to be superlinear but smooth (not abrupt like a piecewise linear characteristic) 19J], and 



second, tunneling makes the junction resistance temperature-dependent Il9| . Since neither 



feature is observed, we can rule out tunneling through pinholes. Refs. [20|, |2l| carried out 
high resolution transmission electron microscopy studies of ferro magnet /organic junctions 
and found them to be abrupt with no evidence of interdiffusion. This further eliminates the 
existence of pinholes in the organic since they would have caused interdiffusion. 

The observed I-V behavior is however very consistent with the transport picture presented 
in refs. At low bias voltages, the Schottky barrier at the injecting contact is thick 

enough to suppress tunneling so that the current is mostly due to thermionic emission. With 
increasing bias, the interface Schottky barrier becomes progressively thinner owing to band 
bending and the tunneling increases. At some threshold bias, the tunneling current exceeds 
that due to thermionic emission. Thereafter, the tunneling injection dominates and with 
increasing bias (decreasing barrier thickness) it increases rapidly. Thus, the current remains 
small up to a threshold bias (~ 2V), at which point cross over from thermionic emission 
to tunneling takes place, and then the current takes off. Since the tunneling probability 
is independent of temperature, the I-V characteristic is virtually temperature- independent. 
Therefore, the observed I-V characteristic is in qualitative agreement with the transport 



picture presented in refs. 



In ref. 



16], we measured the ratio AR/R as a function of electric current through the 



organic at a tem per ature of 1.9 K |16j. In our samples, L « 30 nm. We assume Pi = 0.4 



(cobalt contact) 22], Po = 0.3 (nickel contact) 22], and, «i = a 2 = 1, which is the same 
assumption as in ref. 0, Q. In reality, the spin injection and detection efficiencies are 



never quite 100%, but they can be very high at organic/ferromagnet interfaces - as high as 



85-90% in some cases 



231 ] . Since L S (E,T) is not very sensitive to a>\ or c*2, our assumptions 



regarding these parameters are not critical in any case. 

Equation (6) relates L S (E, T) to AR/R. Hence, measurement of AR/R at various current 
levels allows us to determine the spin diffusion length as a function of current through the 
organic. Knowing the current, we can find the electric field across the organic as follows: We 
apply Ohm's law to find the voltage V across the organic from the relation V = IR, where / 
is the current and R is the measured resistance of the organic. The average electric field in 
the organic is then found from the relation E = V/L. This allows us to determine L S (E, T) 
as a function of E. The electric field in the organic is of course not spatially uniform, but 
the arguments presented here do not require the field to be uniform. 

In Fig. 1, we plot L S (E,T) versus E. The data show that the spin diffusion length 
L S (E,T) monotonically decreases with increasing electric field. That implies that the spin 
relaxation time decreases very rapidly with increasing electric field. To understand this, 
note first that we are operating in the high- field regime where \E\ ^> kT j (eL s (E,T)). In 
our experiment, the average electric field strength \E\ varied between 3.16 kV/cm and 60 
kV/cm, whereas kT j (eL s (E,T)) varied between 220 and 303 V/cm. This puts us in the 
high field regime. In this regime, Equations (1) and (2) yield 

[Ls(E,T)] MghE « »(E,T)E[r s (E,T)) highE 

efjE/T)VE ( P{T)y/E \ ( P(T)VE - A(T) \\ 



kT r- , , , /3(TWE\ / /3(TWE - A(T) \ r 

-jVEh{0, T)tanh f J exp f Pl ;V fcr 1 ] \ [r s (E, T)) high E 



which tells us that the spin relaxation time [r s (E, T)] hi h E will have to drop off 
super-exponentially with the square-root of the average electric field in the organic if 
[L S (E, T)] high E decreases with increasing field. Since that is the behavior of [L S (E, T)] high E 
we observe experimentally, we conclude that the spin relaxation time must have decreased 
very rapidly with increasing average field in the organic. 

The rapid decrease in the spin relaxation time with increasing electric field is once again 
not consistent with hyperfine interactions. To first order, the strength of hyperfine inter- 
action is independent of the electric field. This strength is proportional to the sum of the 



carrier probability densities (squared modulus of the wavefunction) at the nuclear sites 24]. 
An external electric field can skew the carrier wavef unctions in space (as in quantum confined 



Stark effect 



251 ]) and change the interaction strength, but it requires a very high electric field 



to skew the wavefunction appreciably since carriers in organics are quite strongly localized. 
Even in quantum confined Stark effect, where the carriers are relatively delocalized, it takes 
field strengths of several hundreds of kV/cm to change the overlap between electron and 
hole wavefunctions by a few percent. Therefore, we do not expect the hyperfine interaction 
strength to be particularly sensitive to electric field. 

There is a second effect to be considered. Even though the hyperfine interaction strength 
may not be sensitive to electric field, the spin relaxation rate due to this interaction may 
become sensitive because the ensemble averaged spin relaxation rate (1/t s ) is equal to 
J °° def(e) [l/r s (e)] / J °° def(e), where /(e) is the carrier distribution function in energy space 
e. Since an electric field can change /(e), it could influence (1/t s ), but such influence would 
be small because according to the variational principle of transport, a first order change in 
the distribution function induces only a second order change in ensemble averaged transport 
parameters such as (l/r s ) [261 ] . Therefore, [t s (E, T)] high E could not be a strong function of 
electric field if hyperfine interactions were dominant. 

If [t s (E, T)] hi h E did not depend strongly on electric field - as would be the case with 
hyperfine interactions - then that would make the spin diffusion length [L s (E,T)] high E in- 
crease super-exponentially with the square-root of the electric field according to Equation 
(J7j). This rapid increase in [L S {E, T)] high E with electric field is what ref. jfl also predicted if 
hyperfine interaction is the primary spin relaxation mechanism. However, what we find ex- 
perimentally is not a rapid increase but rather a decrease in [L s (E,T)] high E with increasing 
field. This trend alone indicates that hyperfine interaction is most likely not the dominant 
spin relaxation mechanism in Alq 3 . Of course, hyperfine interaction is suppressed in a mag- 
netic field 0,0], and therefore it is possible that the magnetic fields used in the experiments 
of ref. [jij, quenched the hyperfine interaction. Nonetheless, we can say that at magnetic 
field strengths commonly encountered in spintronic applications, hyperfine interaction is not 
likely to be the major spin relaxing mechanism. 

Finally, there are theoretical objections against hyperfine interaction as well. Most or- 
ganics of interest are 7r-conjugated molecules where the delocalized electron states are p z 
orbitals whose nodal planes coincide with the molecular plane. Therefore, contact hyperfine 



7 



interaction should be vanishingly small in them jl|. In some organic semiconductors like 
Alq%, the electron wavef unctions may tend to localize over carbon atoms [281] , whose natural 
isotope 12 C, has no net nuclear spin. Hence, contact hyperfine interaction should typically 
be weak in organics. 

Before we conclude, we point out that the data in Fig. 1 may actually suggest that 
the Elliott- Yafet mechanism jg] is the major spin relaxer. This mechanism has its origin 
in the fact that any spin-orbit interaction makes the eigenspinors of a carrier in the lowest 
unoccupied molecular orbital (LUMO) or highest occupied molecular orbital (HOMO) states 
of a 7r-conjugated molecule like Alq 3 momentum-dependent, so that whenever an electron or 
hole scatters and loses (or gains) momentum, its spin relaxes. Although spin-orbit interaction 
in organics is weak, it may not be so weak as to preclude the Elliott- Yafet mechanism 
altogether. This is the conclusion we reached in ref. 9j as well. 

In order to understand why spin relaxation via the Elliott-Yafet mode could 
make [t 8 {E, T)] hi h E decrease rapidly with increasing electric field (and therefore make 
[L S (E, T)] high E decrease with increasing electric field - consistent with the data of Fig. 1), 
consider the fact that in this mechanism, the spin relaxation rate 1/tJE, T) is roughly 



proportional to the momentum relaxation rate l/T m (E,T) and is given by 29] 

1 A(E) 1 



(8) 



r s (E,T) E g r m (E, T) 

where A(E) is the (electric-field-dependent) spin-orbit interaction strength in the LUMO 
levels for electrons or HOMO levels for holes, and E g is the HOMO-LUMO gap. The 
momentum relaxation rate l/r m (E, T) will increase with electric field E because of enhanced 
scattering, but more importantly, the spin-orbit interaction strength A(E) will also increase 
with electric field. The spin-orbit interaction Hamiltonian is given by the expression 

H so oc (S x pj ■ ff, (9) 

where p is the momentum operator, E is the total electric field that the carrier sees (which 
includes the externally applied field) and a is the Pauli spin matrix. Therefore, A(E) 
should increase with E. This is a well known fact in inorganic semiconductors and one of 
its manifestation is the celebrated Rashba effect [30] . We can expect a similar effect in a 
disordered organic as well. That, taken together with the fact that the momentum relaxation 
rate also increases with E, should make the spin relaxation time t s (E,T) decrease rapidly 
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with increasing E. This then should make L S (E,T) decrease with increasing electric field, 
consistent with our experimental observation and the data in Fig. 1. 

The Elliott- Yafet mechanism is also consistent with the observed temperature depen- 
dence of the spin relaxation time in ref. js]. According to Equation (jSj), the spin relaxation 
rate and the momentum relaxation rate should have the same temperature dependence since 
A(E) and E g are nearly temperature independent. Hence t s (E, T) should exhibit weak tem- 
perature dependence if Coulomb scattering is the dominant momentum relaxing mechanism, 
since this scattering mechanism is elastic and makes the momentum relaxation rate nearly 

n 

temperature-independent. In ref. [9J, the major momentum relaxing mode was Coulomb 
scattering and t s (E,T) expectedly exhibited weak temperature dependence. Therefore, in 
the end, both the temperature and the electric-field dependences of spin relaxation time are 
consistent with the Elliott- Yafet mechanism, but not with hyperfine interactions. 

In conclusion, we have shown that experimental evidence gathered so far and theoretical 
considerations tend to favor the Elliott-Yafet mechanism more than hyperfine interactions 
as the dominant spin relaxation mechanism of carriers in the most widely studied organics. 
Nonetheless, further experiments are required to resolve this issue conclusively. 
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FIG. 1: Spin diffusion length as a function of electric field. The measured spin diffusion length in 
50-nm diameter Alq% nanowires as a function of electric field at a temperature of 1.9 K. 
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